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SUMMARY 


The investigation of the dynamics of global atmospheric disturbances 
is essential for the understanding and control of the earth environment. A 
theoretical model based on the kinetic theory for the perturbation of plasma 
in the magnetosphere is proposed to study the observed disturbances which 
are caused by both natural and artificial sources that generate wave-like 
perturbations propagating around the globe. The proposed model covers the 
wave propagation through a media of transitional (from collisional to 
collisionless) fully ionized magnetoactive plasma. In the present study, we 
have presented a systematic formulation of the problem and the method of 
solution for the transitional model of magnetosphere is discussed. Finally, 
the possible emission of hydromagnetic waves in the magnetosphere during 
the quiet and disturbed time are also discussed. 


ii 



CHAPTER I 


INTRODUCTION 

The dynamical behavior of the magneto spheric disturbances is 
essential for the understanding and control of the earth environment. 

In particular j accurate information about global atmospheric disturbances 
is required for the design and operation of satellites, space shuttle and 
space laboratory, for improved communications and weather prediction, and 
for monitoring and controlling earth resources, all of which are closely 
related to the improvement of our living conditions. 

The observed disturbances are caused by both natural and 
sttificial sources which generate wave— like perturbations that propagate 
through the magnetosphere. Wave— like phenomena caused by natural sources 
refer to disturbances induced by non— man— made sources such as solar flares, 
solar wind, and particle precipitations. The wave— like phenomena caused 
by artificial sources refers to man-made disturbance such as nuclear 
explosions, radio wave propagation, space vehicle, rocket launches, etc. 

The purpose of the present research is to seek a systematic investigation 
of the wave— like disturbances generated by either natural sources or 
artificial sources which originate in the magnetosphere. In particular, 
it is of considerable interest to study how waves are generated and behave 
and how the wave-like disturbances are dissipated. 

In general, propagation of waves or wave-like disturbances is 
modified by transport phenomena due to Coulomb collisions. Collisional 
effect can vary from region to region in the medium (i.e., fully ionized 
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plasma) of interest. In particular, waves of a given, period may see one 
region of a medium as collisionless (in the sense that the wave period is 
short compared with Coulomb collision time) and another region of the 
medium as collision-dominated. When a wave propagates upward in the iono- 
sphere, the Coulomb collision frequency decreases. Thus, for the wave or 
wave- like disturbances of a given period, transition from collision to 
collisionless behavior may take place. 

Transport coefficients of Coulomb collisions are greatly modified 
by magnetic fields. Thus, when the waves or wave-like disturbances propa- 
gate in a magneto-active plasma, the following modification should be con- 
sidered (Braginskii, [1]); 

(1) S2./V, ( Syrofrequency \ ^ magnetic 

^electron collision frequency / 

field. 

(a) For the case of waves propagating transverse to magnetic 
field (B), the magnetic field strongly modifies thermal 
conduction and viscosity. 

(b) For the case of waves propagating parallel to the 
magnetic field, there is no effect on transport 
coefficients . 

(2) ^ » weak magnetic field 

In this case, the magnetic field does not modify the transport 
coefficients for waves propagating either transverse or 
parallel to the magnetic field. 
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(3) ~ 1 

In this case, the modification of transport coefficients 
is very complicated. It needs a full solution of the 
kinetic equation with magnetic field effects. 

In this investigation, we shall present a systematic transitional 
model for the studying of propagation of waves or wave-like disturbances in 
the magnetosphere. The model is based on the kinetic theory, and the method 
of solution is followed by the work given by Hung and Barnes [2], [3], and [4]. 
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CHAPTER II 


FORMULATION OF THE PROBLEM 
II“1 Basic Parameter 

As we have discussed In Chapter I, the behavior of wave-like 
disturbance propagation is modified by transport phenomena due to Coulomb 
collisions, and the collision effects are determined by the ratio of the 
collision-frequency to the wave-frequency (v/co). Suppose that an observer 
fixes his system of inertia coordinates with a wave of a given wave- 
frequency oj, and observes another moving coordinate system fixed with a 
coordinate system with the collision frequency V. The observer in the 
coordinate system co reaches the conclusions concerning the particles motion 
as follows: 

(1) when CO « V, particle motion is a continuous motion, 

(2) when co ~ V, particle motion is in transition from continuum 
to discrete motion, 

(3) when co » V, particle motion is resembling as discrete motion. 

From the above mentioned facts, we may classify case (1) as being 

collision-dominated; case (2) as being transitional from collisional to 
collisionless; and case (3) as being collision-free. Thus, for studying 
the dynamics of wave-like disturbances of a given period, it is obvious 
that the selection of appropriate physical models depends on the ratio of 
collision-frequency to wave-frequency (v/w). Therefore, the basic parameter 
for studying such a problem will be v/w. 
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In the magnetosphere, the collision time for electron and ion can 


be calculated from 


3.5 X 10^ ^e ^ 

-r? (sec) 


X 


m 


( 1 - 1 ) 


^ 3.0 X 10^ /™i 

^ ^ * Zra^ 


1/2 m 3^2 


T. 

1 


X 


nn 


(sec) 


(1-2) 


with mp and Z being the proton mass and atomic number of ion, n and X are 
defined as 


n = n^ = Zn^ (1-3) 

and X is called Coulomb logarithm, 

{ 23.4 ~ 1.15 iln n + 3.45 Te < 30 ev. 

(1-4) 

25.3 - 1.15 V n + 2.3 Tg for > 50 ev. 

As soon as we have decided the wave-frequency of our interest, the 
v/o) can be calculated, thus the mathematical model for a particular problem 
can be followed. The details of this mathematical model are presented in 
the next section. 

II- 2 Governing Equations 

During the present study, we have been interested in the propagation 
of ULF, ELF, and VLF waves in the magnetosphere. Their frequency ranges are 
3 X 10'^ Hz - 3 Hz, 3 Hz - 3 K Hz and 3 K Hz - 30 K Hz, respectively. 

Most ULF waves can be identified as Pc 1, Pc 3, Pc 4, Pc 5 (quiet time 
transverse); Pc 4, 5 (quiet time compressional) ; Pi 1 and Pi 2 (substorm 
compressional) ; Pc 1 (storm transverse); Pc 4, 5 (storm compressional). 

For this case, we can idealize the problem by treating the electrons as a 
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fluid (collision-dominated case) and assume that the ions are adequately 

described by the Boltzmann equation that neglects ion- ion, but not ion- 

electron collision. The validity of this idealization has been discussed 

by Hung, Wu and Smith [5]. Also, we have considered that the plasma and 

magnetic fields are uniform, on the average, throughout an effectively 

infinite volume, and there is no average electric current, and the 

average electron pressure tensor is isotropic, but the ion pressure 

tensor is not. Let n, v, T, p. q and j denote number density, velocity, 

temperature, pressure, heat flux, and viscous stress tensor, respectively. 

The superscripts and subscripts i and e denote electrons and ions. The 

other symbols are; E, the electric field; B, the magnetic field; e, the 

ion charge; i the ion or electron mass; and c, the speed of light. Then, 
(e) 

the fluid equations for the electrons can be written as follows; [1] 


-^ + V ■ (n^ v^) = 0 

at 


( 2 - 1 ) 


'•Xe - 


dt 


= - V - V • TT^ - ene (E + v® x B) + 


( 2 - 2 ) 


3 die 

7 >^6 TT + Y-f 


(2-3) 


The kinetic equation for the ion velocity distribution f^ is 


3£i e 1 

i + V V f. + |- (E V X B) • 


9t 


m 


^3 Te 3fi 

• (vf . + — — ) - 

^i^e ~ ^ ®i ^ 


R. 

~x 


1 “1 


9v 


(2-4) 



where 


P 

e 


n T 
e e 


( 1 ) 

P. II 
1 


(J-) 


= n.T. 
1 1 


(2-5a) 


^ ^ -i. u 

dt “ IF ~e ' ^ ’ (2-5b) 

and T is the electron collision time, and R = -R. is the collisional 
e "^e ~x 

momentum transfer from ions to electrons, R is composed of a frictional 
force R and a thermal force R , in which 

Vi ^ X 


R = R + . 

~u ~T 


( 2 - 6 ) 


The electron thermal flux is composed of analogous parts, 

_ u T , 

9e ^ 9e’ shown from the velocity moments. Finally, of course, 

the electromagnetic fields must satisfy Maxwell's equation 


V • E = 

V X B = 


4TTe(n.-n ) 
1 e'^ 



1 ^5 

c 3t 


V X 


E 


9B 
c 3t 


= 0 


V • B = 0 


(2-7) 


These equations are in Gaussian units. 



CHAPTER III 


METHOD OF SOLUTION 

III-l General Procedure 

In this section, we shall outline a general procedure to obtain a 
solution for this set of governing equations, (2-1) through (2-7). Due to 
the complexity of mathematics in nature, it is formidable to obtain an 
analytical solution from this system. However, the physical characteristics 
of this present problem process a wide range of variety. Before we obtain 
a complete solution, it is still necessary to attract some important 
physics from this complex system, such as the hydromagnetic Instabilities, 
dissipation rate, etc. Therefore, we shall proceed to examine those meaning- 
ful parameters. 

As discussed in Ref. (5), the time scale of interest in the present 

theory of transitional model is that electrons are in transition from 

collisional to collisionless and ions are in collisionless. Assuming that 

ion and electron temperatures are in the same order of magnitude, it was 

found that the ions become collision-dominated when the wave period 

1/(a) ^ (mi/mg) (where (a, m, and denote wave frequency, mass, and 

ion collision time, respectively); and electrons become collision-free when 

1/oj ^ T^. In the present theoretical model, the electron equation 

breaks down when l/o) <. Tg (where electron-electron collisions become 

insignificant); and the ion equations are invalid when l/u) > (m /m T 

— 1 e i 

(where ion-ion collisions play a significant role). Therefore, the present 
theoretical model is applicable only for the time scale (or wave period) 
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As it stands, Eqs. (2-1) - (2-3) show the governing equations 
of electrons, and Eq. (2-4) is the governing equation of ions* Eq. (2-2) 
includes the collisional momentum transfer from ions to electrons, and 
we have Ignored collisional energy transfer ions to electrons in Eq. (2-3). 
This is because the collisional ion-electron energy exchange, whose 
characteristic time is on the order of or greater than (ra^/m^)^^^ 
which is greater than the time scale of present interest. 

has been neglected. The right-hand -side of Eq. (2-4) shows the ion-electron 
collision term which is of the same form as the Fokker-Planck collisional 
term that describes random motion of particles in a moving medium with 
temperature . The first term of the collision terms describes the 
collisional energy transfer from electrons to ions, and the second term 
implies the collisional momentum transfer from electrons to ions. Again 
we are going to neglect collisional energy transfer term because the time 
scale of the collisional energy transfer which affects the evolution of 
distribution function is on the order of or greater than Tg(m^/m 0 ) or 
Ti(mi/mg)^/^ that is long compared with the timescale of current problem 
of interest. 

In the present analysis, we are concerned with waves whose 
circular frequency w is small compared with the ion gyrofrequency , and whose 
wavelengths are long compared with the mean ion Larmor radius. Under such 
circumstances, the momentum transfer due to collisions can be represented by 

= 0.71 V|| Te (3-2) 

*Here the unit of T is the erg. 
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where the subscripts ||, J_ refer to the magnetic field direction e^ = B/ j 
Similarly, the electron heat flux is 

9e = -Ke!|l «-3) 

where Kg is the coefficient of electron heat conductivity. Furthermore, 
the stress tensor after Braginskii [1] is 


-n-e -qC . 

where the rate of strain tensor, w^, is 


(3~4) 


w 


8v^ 
e _ a 


a3 


+ e. - 1 6 


8xc 


3x, 


9v® 
a 

3 '"“6 ^ 


(3-5) 


a 


and n is the tensor coefficient of electron viscosity, which is the 

J 

function of and and subscripts a and 3 represent the coordinates. 
Under the present condition fig Tg » 1, the stress tensor has the 
following form in a coordinate system with z-axis parallel to the magnetic 
field 


ZZ VO 


8v® 

z 

3, 


9z 

3 

V 9x 3y 3z /_ 


(3-6) 


n = II = HI 
XX yy VO 


3v® 

x 


8v^ 


3x 


+ 


!Zz 2 _ 

3y - 3 \ 3x + 3y 3 


3v® 

z 


(3-7) 


TT^ ~ ~ IT® 

xy yx xz 


e e 


zx 


yz 


lT®y ~ (l/fig/T^) 
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(3-8) 


where the zeroth-order coefficient of electron viscosity is 


n = 0.73 T 

VO e e e 

For the convenience, we shall seek solutions for electrons and 
ions separately. The details for these manipulations will be presented 
as follows: 

(i) Electron Dynamics 

In order to solve the fluid-like electron equations, we assume 


T 


<T> + 6t (x, t) 

n 


<n> + 6n (x, t) 

P 


<P> + 6P (x, t) 

B 


<B> + 6 b (x, t) 

E 


5 E (x, t) 

e 

V 


6 (x, t) 



— 


where < > denotes, ensemble averaging. We consider the limit of small 
amplitudes fluctuations, j5n / <n> | << 1, etc. If the fluctuations are 
sinusoidal, i.e., proportional to exp [i (k • x - u?t)] where i = 

Eqs. (2-1 through 2-3) become after linearization and neglecting 
terms of order m^/m^ 


(a) 6n^ 


<n^> 


(k I 6v^ 
1 X 




(3-10) 
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1 6 Pg k I + Pvo (I “ f ’^1 ‘'I I ® 


(3-11) 


+ m <n > ^ 6v =. 0, 
e e e y 


e6E - in ^ 6v = 0 , 
y e e X 


(3-12) 


e 1 


16 P^ k^l + 2tVo (I k| I «v® - j kj^ k| I 6v®) + e <n^> «E„ 

(3-13) 


+ 0.71 i <n > k 6T = 0, 
^ e 


~ " (f **■ ^ ^| | ^e) • (3-14) 


Here, without loss of generality we have assumed k = (k^, 0, ^||)‘ Combining 
Eqs. (3-10) and (3-14), we have the following relation for pressure and 
number density 


«Pe = 6 


where 


Y = 1 + 

e 3 4- 2 i n 


(3-15) 


(3-16) 


TC 


Htc 


ki I K 


0) <n^> 


(3-17) 


To Investigate the physical significance of Eq. (3-16), let us 
examine the parameter T in the following way: 


p = ^ while K ^ ^ 
3 e 


r := 1 while K 

e 


(3-18a) 

(3-18b) 
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This means that in the two extreme cases of zero and infinite thermal 
conduction, the electrons fluctuate adiabatically and isothermally, 
respectively. 

The velocity fluctuation of electrons can be easily found in terms 
of 6 E by using Eqs. (3~10) and (3-15) in Eqs. (3-11) through (3-13) 


'6v®' 

" 

X 


6v® 

c 

y 

""<B> 

6v® 


z 

•• — 



0 


M 


xy 




yx 






M 


yz 


m: 


zy zz 

where the components of the mobility tensor are 


" 

'SE^' 


6Ey 

J 



(3-19) 


== = 1 
XV yx ’ 




kf a^ C 
- 1 ^ ^ 

_ 


1 + 


2(3r^ + i2c 




2i 0) 




zz 


2i.2 


(1.71 r -0.71)a^k 
e e 






zy 




4i^1 

1 + 

L “ J 






yz k 


a. 


where a = 3 (1.71 - 0.71) - ^iC^> and 


con 


r (oi) = — = 0.73 COT 

^ <p > 


(3-20) 
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is in general a complex function of co, and = (2 is the 

electron thermal speed. From Eqs. (3-20), (3-18a, b) and Braginskii^s 
expression for K^, it follows that 


where 3 = 8tt <n> 

G 

Alfven velocity. 


n 


TC 


m. 3 

2.16 ~ -I C 

m u V 


(3-21) 


e 

<T >/<B>^, u = (A)/|k,j|c , and C = <B>/(47itn. <n>)^'^^ 
^ II A A X 

In deriving Eq. (3-19) we have exploited the fact that 


< < 1. 


(ii) Ion Dynamics 

Next, we consider the ion kinetic equation, with the object of 

finding an expression analogous to Eq. (3-19), for ions. The first terra 

on the right hand side of Eq. (2-4) affects the evolution of the velocity 

distribution on the time scale t m./m - t . (m./m )*^^ which is long 

e X e X 1 e 

compared with time scales of importance for the wave, and may therefore 
be neglected. As we pointed out earlier, the frictional force is also 
negligible. Hence, from Eq. (3-2), the right hand side of Eq. (2-4) is 
just 


1 • 

i 

m. n. ~1 3v 

XX 


0.71 

m. 

1 


(V 


e 



(3-22) 


Linearizing Equation (2-4), we obtain 

0<f,> 

(v X <B>)- = 0, (3-23) 


whose solution is 

<f , > = f (v I I , 
X O I I 

where f^ is arbitrary, and 



(3-24) 
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1-^ + (v “ V) + — “ (v X <B>) • 6f . 

ydt - m^C - ~ 9 vJ i 

3<f.> 

"lyjr (3-25) 


Finally, the Faraday's law becomes 


= - — ( 6 e + — V X 6 B) 
~ c - 


9<f > 

. 0.71 , 

_ 1 k|| oT • — 5 — — 

m. - e dVii 


~ 6 B = k X 6 E. (3-26) 

c ~ ~ 

Now, we have reduced the governing equation to a workable form (i.e,, 

Eq. (3-10) through (3-26)). In the following section, we shall describe 
the procedures to obtain a dispersion relation from this set of equations. 
III-2 Derivation of Dispersion Relations 

Eq. (3-10) through (3-26) permits us to write 6 b and as linear 
combinations of the fluctuating electric field components. The Eq. (3-25) 
can be solved by standard techniques [ ] for giving an ion mobility tensor 
analogous to the electron mobility tensor of Eq. (3—19). However, it is 
simpler to exploit the fact that Eq. (3-25) is equivalent to the 
linearized Vlasov equation with 6 E replaced by 6 e + i (0.71/e) kjj 6 T^ e^. 
Hence, if <St^ is expressed in terms of 6 e, we may easily find the ion 
mobility tensor from Vlasov mobility tensor. 

From Equations (2-14) and (2-18) we have 


(r - 1) <T > 

(ST = — k “ (M ) • 6 e 

e <B> ~ ~e 


By using the expression from (A-23) we can show 
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where 


and 


Thus, 


k • (M ) • (5E 


= (k^, 0, k| ij 


-1 


k? 


1 e 

1 — o ^ 
lx)U 


1. 

1 kii ""2 


-j— — D- i ’ T t ' 2 ~ D / 

ki I 3 a k I 4 

11 e I 


0, -i DC^k|, 



, 3T + i 2^ 

_ 1 . e 

D = _ + 

12 a 


°2 = 


jr +12? 

e ^ 

a 


D = 1 + ^ 

3 a 


iH, 

D, — 1 + 

4 a 


D = 


3(1.71 r - 0.71) - i4C 
e V 


we have 


3i(r -1) 

6t = 

e a 


k 1 a' 


1 

k|t ~ (M 
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Making use of Equations (3-9), (3-24) and (3-27), we rewrite 


Equation (3-25) 


[-iw + ik*v+J^^ (vx e^) • -^] 6f ^ 


1 I 

~ (6e +-^ v X 6 b) * — A I (i 

_ 


' \ 3<f.> 

«E - —t; ^ ?. «E._ -5-i- 


(liQ^ -V y J 9v j 


L -^\ 3 - 28 ) 

where A = 3(r^-l)/a. The coordinate system is chosen such that k = (kj^, 0, kjj) 
and V = (vj^ cos\|;, vj^ sin4^, V| |). Then Eq. (3-28) is an ordinary differential 
equation in if;. 


1 1 


dT“ ^ ' ll I 


(k| I Vi I - iu) + k| V| cos Ip) 6f . 




= - <f .> 
1 


kiki la' 


6e / Tlr^ 

77” I V I cosif; — T + V I sinif; \ 1 + A t, 1 — ^ (3-29) 

i <Tt> 1 \ <tS 


+ VI I (1 -A) 


<Tm> 


The formal solution of the ion kinetic equation (3-29) can be solved in 


terms of Greenes function, i.e,. 


6f - 

^^i n 


— f V| cosif; r + V| sinib [l 

ij |1 <Tt> i V 


l ll e \ 6F, 

+ A ^ ^ 

“^e 7 <t]-> 


(3-30) 


+ V„ (1-A) — ^ <f > G (r.V) 

I! <T,,> 


where the Greenes function satisfies the particle path along the non- 
perturbed orbit, namely 

G(i|;; if;’) == exp [-w (if^ - if;’) + kj^ vj^ (sinif; - simf;') 


+ k| I V. I (if> - 


(3-31) 
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By using the relations 


1 

cos 

sin 


exp (-i A sini|^) = 

n=-oo 


Equations (3-30) and (3-31) give 


Jn (A) 

- Jn (A) 
i Jn* (A) 


exp(-in^Jj) (3-32) 




i exp (i (kj^ - n i/; ] } 


n= 


kj I VI I - w + nf^^ 


2 




2!^^. 

n J ~ dE +i ^ 
n <B> X 


ki kiia' 


— - — V I J 


v,J 1 + A 


1 




C (3-33) 


2 ^. 


6E + 


<B> y 


(all)" 


\ (1 - 


<B> z 


<f .> 

1 


where J =J (kiVi/f^.) and J are the Bessel functions of the first kind 
n n J_ ^ X n 

and order n, and its first derivative. 

Maxwell’s equation gives the relation between electric current 
density j and dielectric tensor K, i.e.. 


j = ^ (1 - J) • 


(3-34) 


where the current density and dielectric tensor are 


j = E e 

j 


. 1 V f. 

J J - 1 


dv 


j = e and i 


(3-35) 


and 


K = — ^ e <n> (M. - M ) + 1 
~ 0) <B> ssX !«e w 


(3-36) 
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respectively. Here is mobility tensor which satisfies the relation 


6v, = -rrr M. • 6 e 

-J <B> - 


(3-37) 


where 6v. is defined 
^3 


6v, = r V f.dv/ jf.dv 

J - J - j J - 


(3-38) 


It is clear that in order to obtain dielectric tensor, we have 
to calculate the current density from distribution function. In the 
present stage, we have to calculate ion current density which follows 
the formula 


e <n> 5v. 


j.=ejvf. dv = 

~i j ~ 1 ~ 

Substituting Eqs. (3-24) and (3-33) we have 


(3-39) 


6v 


6v^ 

y 


6v"" 

z 


=-i 


00 oo oo 27T 


vj^cosij; 

Vj^sinijj 


(ki|Vjj-(jj+n^^^) 


1 


■irn 


• exp 




ktki la" 


1 + A 


i 




c I + <Se 

v/ z 


nJ + 6E 

”1 > 


2^. 

7 r ,2 I 

(an) ' 


2a. 

— i— V| j’ • 

, i 2 i 

(ap 

2 

1 


exp 


(aj^) 


i 


(at. ) 


(3-40) 


By using the relations 
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2tt 


A. 

2tt 


27T 




1_ 

27T 


0 

2tt 

f 

0 

2tt 


± 1 (z sin X - nx) , a. , / \ 
e dx=±J (z) 

n 


+ 1 (z sin X - nx) , . n ^ . 

cos X e dx = ±“ J (z) 

z n ' 


sin X e 


± i (z sin X - nx) 


dx = — i j’ (z) 
n ' 


(3-41) 


Equation (3-40) becomes 


“1 



nO, . 


CSO 00 ( 


J 

^3y 

i r ^ \ Tl 

r 

dV|| 

k|V 1 n 

~ T 

-1J„ 

6ji 

_ 

0 —00 

J 

n 


2U. 

^ T ^ 

nJ , 


, 1.2 n* > 1.2 n 

(aj_) (aj_) 


\ 2!;^. 

1 + A C , — ?^v,,J(l-A) 

ojQ ■ V / ’ /lx I 1 n 
(ail) 
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(3-42) 


Under the hydromagnetic assumption, the wave frequency w and bounce 
frequency k| j V| j of typical particles (ions in this case) are well below 
the ion cyclotron harmonic frequencies, i.e.. 
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Let us expand the summation of ion cyclotron resonance denominators in 
Eq. (3-42) into the following forms 
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By using the following relations 
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we have the following expansions 
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Making use of Eq, (3-46), Eq, (3-42) leads to 
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Let us introduce the plasma dispersion function 2(J) which is defined 
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In our case, the following relations in terms of plasma dispersion functions 
are most useful 
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Then we can include contour of integration into the plasma dispersion 
function which makes Eq. (3-47) become 
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where 
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By using the following integrations 
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where the mobility tensor with anisotropic ions are 
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where A = 2,13 (F - l)/ot, = e <B>/m. c, i\ “ (2 <T^>/m.) 
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= w/ (|kj j ]a I j ) , and = Z(y^) and are the plasma dispersion 
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function and its first derivative (Fried and Conte 1961). In deriving 
Eq. (3-53) we have neglected the terms proportional to A because 




(equation continued) 
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The requirement that Eq, (3-55) be consistent with Maxwell 
equations gives the dispersion relation 
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Substituting Eq. (3-57) in Eq. (3-58) gives 

2 


ck| 


XX \ to 




k|c^ 


zz 


o 

-*• kT 

2 / yz zy 


= 0 (3 


Eq, (3-59) is given at the lowest order in w/f2,. The first factor in 


Eq. (3-59) is the usual Alfven wave dispersion relation 
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(3-60) 


and the second factor, which gives the magneto-acoustic dispersion relation, 
is 
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since |k^z I ^ ^ |k c/w| ^ . 

By using Eq. (3-57), we can rearrange Eq. (3-60) by straightforward 
calculation to give 
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and u = 00 / ( I k I I I C^) . 


(3-69) 


We notice that if dispersion relation (3-62) is satisfied for a 
wave (o), k) , it is also satisfied for a wave propagating in the opposite 
direction, i.e., for the waves (ia), -k) and (~w*, k) . This is so first, 
because (3-62) is invariant to the transformation (a), k) ^ (w, -k) and 
secondly, because the transformation ( 03 , k) k) causes 

C - Z * 9 r r* , Z. Z’. *, (u*)^, so that Eq. (3-62) is 

transformed into complex conjugate. Hence, in particular, we may always 
choose Re > 0 and R e ^ > 0. 

Note the superscript * Indicates the complex conjugate. 
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CHAPTER IV 


MAGNETOSPHERIC HYDROMAGNETIC INSTABILITIES 
IV- 1 Introduction 

Observations from Explorer 12, 14, and 15 indicate a large amount 
of proton precipitation in the magnetosphere with energies 100 KeV to 
10 MeV [7]. In the meanwhile, data from Explorer 12 confirms a large 
flux of trapped protons in the energy range 100 KeV to 4.5 MeV with 
peak intensity near geomagnetic shell L - 3.5 [8]. By using the 
parameter at L = 3.5 (intensity of geomagnetic field ~ 10“^ gauss, proton 
number density ~20 cm”"^, energy of protons ~200 KeV), the ratio of 
plasma pressure to magnetic pressure 3 becomes to the order of unity. 
Therefore, it is interesting to Investigate the emission of hydro- 
magnetic waves due to possible plasma instabilities, the dissipation 
of hydromagne tic waves and how it may correlate to the modification 
of heating mechanism as 3 approaches unity in the magnetosphere. 

Experimentally, several authors, [9], [10], have observed that 
hydromagnet ic waves in the frequency range of Pc-I pulsations can 
be generated by proton beams traveling faster than the Alfven velocity, 
if the beam has a certain anisotropy in pitch angle. This implies that 


a . 
1 


2T^ /m . 
i 1 


B^/47Tm.n 

1 




(4-la) 


where a. is the proton thermal velocity; C . , the Alfven velocity; and 
n, the total number density, Jacobs [11] estimates the hydromagnetic 
wave with experimentally observed characteristics can be generated at 
geomagnetic shell L ~ 4.0 by protons with energies 200 to 500 KeV. As 
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we have mentioned earlier, the observations based on Explorer 12, 14 
and 15 confirm the estimation. The hydromagnetic waves can also be 
generated at L 5-7 by protons with lower energies (several tens of 
KeV) , At L = 3.5 and L = 7.0 the proton gyrof requency , is 96 

and 8.6 Hz, respectively. As has been discussed by Hunt, Wu and 
Smith [5], the time scale of interest for the study of hydromagnetic 
waves is chosen to be 0.2 to 10 seconds. By using this time scale, 
the ratio of hydromagnetic wave frequency to the proton gyrof requency, 
is much less than unity at L < 5, and on the order of or greater 
than unity at L > 5. In the present study, we limit ourselves to 
investigating the hydromagnetic instabilities which could be relevant 
to the emission of hydromagnetic waves in the magnetosphere at a 
distance L < 5, while ion cyclotron instability could be responsible 
for the emission of waves in the periods of 0.2 to 10 seconds at 
a distance L > 5. The model of ion cyclotron instability in the 
magnetosphere has been investigated by Cornwall [9] , Feygin and 
Yakimenkd [12] and Gendrin, et al, [13]. 

As has been mentioned in Ref. [5], the propagation of hydromagnetic 
waves is modified by transport phenomena due to Coulomb collisions, 
Collisional effects can vary from region to region in the plasma of 
interest. At L = 3.5 the electron collision time is on the order of a 
second, while the ion collision time is on the order of a minute in 
the magnetosphere. Thus, for wave periods of 0.2 to 10 seconds, magneto- 
spheric electrons are in a transitional regime between collisional and 
collisionless conditions while ions are in a collisionless regime. 
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IV“2 Criteria for Mirror and Fire-Hose Instabilities 

Most of the known plasma instabilities, which are thermal over- 
stability, are a kind in which (jj is the complex rather than pure 
imaginary [lA], The non-occurrence of overstable solutions has been 
studied by Vedenov and Segdeev [15], Chandrasekhar, et al. [16], and 
Barnes [17] for the case of collisionless plasma* The significance 
of the nonexistance of overstable solution is simply that the con- 
dition of marginal stability for the plasma is given by the dispersion 
relation with w=0. This conclusion remains valid for the case of 
collisional electrons and collisionless ions. 

Let us rewrite dispersion relation, Eq. (3-62), in the following 

form: 


D (u, 0, = j^l+l (6| - 3||) - U^J cos^e 

(4-1) 

-S(u, WT^) sin^6 . 

Assume that w and 0, or u and 0 are an overstable solution of Eq. (4-1). 
Then 

Im [ D (u , 0 , WT ) ] 
e 



(4-2) 
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with = z(y^), and being the plasma dispersion function and its 
first derivative [18] respectively. The y^ is given by 
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Hence, if sin 0 0, I id (D) cannot vanish unless U = 0, so that Eq. (3-62) 
has no overstable solution unless sin 0 = 0. If sin 0 = 0, the dispersion 
relation reduces to Eq. (3-60). 
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The nonexistence of overstable solutions recovers the familiar 


fire-~hose and mirror instabilities for collisional electrons (isotropic 
pressure tensor) and collisionless ions (anisotropic pressure tensor). 
The criteria of the instabilities in the present case can be expressed 
as follows: 


i 
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S|! 


< -1 


(4-4) 




< 0 (4-5) 


Equation (4-4) is the Alfven fire-hose instability to which Equation (3-60) 
is relevant, and Equation (4-5) is the magneto-acoustic fire-hose 
instability to which Equation (3-62) is relevant. In both cases, there is 
an angle 0 q at which the stability is marginal. For the magneto-acoustic 
fire-hose Instability unstable waves propagate at angles 0 such that 


0 <. 0 < 0Q or ir - 0Q < 0 <. 7T, (4-6) 

and for the mirror instability unstable waves propagate at angles 0 such 
that 

^<0<7r - 0Q (4-7) 

The equation for Oq is just D (0, Oq, 0) = 0, which may be rewritten as 
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(4-8) 


This angle 8^ is real only if one of the instability criteria (4-4) and 
(4-5) holds, 

IV- 3 Instabilities Under the Magnetospheric Conditions 

Numerical analysis of the dispersion relation (3-62) , which is 
relevant to the propagation of hydromagnetic waves, has been analyzed 
under ionospheric conditions by Hung, Wu and Smith [5], In the present 
study, we limit ourselves to the investigation of hydromagnetic 
Instabilities which could be relevant to the possible magnetospheric 
disturbances. 

For the plasma with anisotropic ions, ^ we have 

shown that the mirror and/or firehose instabilities might occur based 
on the instability criteria (4.4) and (4,5). The former is the Alfven 
firehose instability in which the triggering mechanism is 

> ej- , (4-7) 

and the latter is the magneto-acoustic firehose instability in which the 
setting up mechanism is 

> e||, (4-8)* 

*In order to satisfy the criteria (4-5) , the second term of the equation shall 
be negative. Since the square— bracket is always positive, we have to make semi- 
circular bracket to be negative. 
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provided that 3^ is on the order of unity for both cases. In other words, 
the instability criteria shown in (4-4) and (4-5) indicates that for the 
case with cold electrons and hot ions, which is typical in the magneto- 
sphere, the instabilities can be triggered for anisotropic particle 

i i i i 

distributions when either 3j | > or ^ ^| | that a strong 

diamagnetic effect for ions (indicating 3^ 1) exists. 

The energetic particle velocity distribution in the magnetosphere 
is, in general, anisotropic. This anisotropy may be caused, for example, 
by charged particle trapped within the magnetic mirrors. The mean energy 
of movement of the charged particles which remain in the trap across the 
lines of force of the geomagnetic field should be higher than the mean 
energy of particle movement along these lines. The experimental data [19] 
Indicates that such an anisotropy of proton velocity distribution in the 
magnetosphere does exist. The observations made by Explorer 26 [20] also 
show that ~ ^ ^ 1 in the magnetosphere. These observations 

show that the magnetoacoustic firehose or mirror instability is one of the 
candidates for triggering hydromagnetic waves. 

On the other hand, the Alfven firehose instability shown in (4-4) 
also is a strong candidate for exciting hydromagnetic waves. This is 
because the dynamics of energetic ions overwhelms the cold electrons, 
and the ions are so nearly collision-free that the ion dynamics are mainly 
governed by Landau damping in which the wave-particle resonant interaction 
changes only the logitudinal component of the particle energy while 
leaving the transverse component unchanged [21]. In other words, 3| | 
is always greater than 3j^ under ion Landau damping conditions when the 
interaction starts with 3| j ~ 3j^. This means r.hat the Alfven firehose 
instability must be considered in the present study. 
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The triggering of hydromagnetic waves in the magnetosphere depends 

on two conditions, i.e., (1) strong diamagnetic effect and (2) 

i i 1 i 

anisotropic properties, either ^ | or 8| | > The latter con- 

dition can be easily satisfied through either the magnetic mirror effect 
or ion Landau damping. This means that the key point for triggering 
hydromagnetic waves depends strictly on satisfying the strong diamagnetic 
effect condition, i.e., B^ shall be on the order of unity. Physically, 
this means that in order to trigger hydromagnetic waves, the pressure 
of the energetic particles in the magnetosphere shall be at least on the 
order of the local geomagnetic pressure. 

It is interesting to compare our proposed triggering mechanism 
with recent hydromagnetic waves observations. Triotskaya and Gul'elmi [22] 
and Jacobs [11] indicate that the propagation of hydromagnetic waves 
become active 1-2 hours before, and 4-7 days after a geomagnetic storm. 

In the meanwhile, more than 50% of hydromagnetic waves are observed when 
the geomagnetic index Kp is less than 2. To correlate these observation 
facts in terms of our model, let us recall the inter-relation between 
geomagnetic storms and solar wind disturbances. It is known that the 
magnetosphere transforms the energy carried away by the disturbed solar 
plasma from the sun into the energy of geomagnetic storms [23]. When 
the enhanced solar wind interacts with the magnetosphere, the plasma 
pressure of the energetic particles in the magnetosphere increases 
dramatically in order to balance the impact from the disturbed solar wind, 
and then the magnetosphere converts the momentum and energy from the 
enhanced solar wind into geomagnetic storms. In other words, the plasma 
pressure is on the order of or greater than the geomagnetic pressure when 
the enhanced solar wind interacts with the magnetosphere, and the plasma 
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pressure becomes smaller than the geomagnetic pressure when the geo- 
magnetic storm develops. The most significant reasons to explain why 
3 is decreasing in the period of geomagnetic storms are as follows: 

(1) the amplitude of geomagnetic disturbances Increase drastically, 
and (2) the observations show that the number density of plasma particles 
decrease dramatically in the magnetosphere [24], [25], [26]. After the 
geomagnetic storms, observations also indicate that a recovery from a 
storm time depletion of magneto spheric concentration takes 4-7 days 
[26]. In the meanwhile, during these recovering periods, the storm energy 
is gradually transformed into Landau damping which heats up the ion 
particles. This dissipation of storm energy into particle energy and the 
balancing of the plasma pressure and the magnetic pressure could take 
several days. This explains why the observed activity of hydromagnetic 
waves increase 1-2 hours before, and 4-7 days after a geomagnetic storm. 
Furthermore, the index is a measure of geomagnetic conditions. This 
means that the geomagnetic pressure is greater than the plasma pressure 
(3 < 1) when the index is high; thus, most of the hydromagnetic waves 

are observed when the K index is low. 

P 

The attenuation rate of hydromagnetic waves propagating in the 
magnetosphere and ionosphere has been calculated numerically based on 
the dispersion relation (3-62) . To match the conditions in magnetosphere 
and ionosphere. Bis chosen from the order of unity to the order of 10~** . 
Thus, the dissipation rate of hydromagnetic waves in ionosphere is 10”^ 
(with = 10”**), and jumps to 10"^ (with = 1) in the magnetosphere 
(see Figure 1). During the periods of higher solar activity with large 
amounts of precipitated energetic particles, lii the ionosphere could 

a 

increase to 10” which makes the dissipation rate of hydromagnetic waves 
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become 10“**. In such a case, the ionospheric heating rate through the 
damping of waves is on the order of 10”’erg-cm”^“sec”^ which is less 
than one percent of the heating rate due to the incident flux of extreme 
ultraviolet (EUV) solar radiation, since the amplitude of the hydro- 
magnetic waves in the ionosphere is only on the order of 1 gamma. Con- 
sequently, we have to conclude that the dissipation of hydromagnet ic 
waves in insufficient to modify the heating of ionosphere even during 
an active solar cycle with large amounts of precipitated energetic parti- 
cles. These results agree with recent observations made by Sorenson [27]. 

On the other hand, the plasma pressure becomes on the order of 
magnetic pressure in the magnetosphere,' This makes Landau damping of 
hydromagnetic waves increase drastically. Our result shows that the 
dissipation of hydromagnetic waves with amplitudes of 10 gamma, which 
contributes to the magnet ospherlc heating rate through the wave damping, 
to be on the order of 10“** to 10”® erg-cm“^-sec"^ which is on the same 
order as the heating due to EUV solar radiation. Hence, we may con- 
clude that the dissipation of hydromagnetic waves could contribute to 
the magnetospheric heating but not to the ionospheric heating. 

In conclusion, we propose that the Alfven firehose and magneto- 
acoustic firehose instabilities could be relevant to the emission of 
hydromagnetic waves in the magnetosphere at locations at a distance of 
L < 5, Our justification is that for a distance of L < 5 the wave 
frequency of hydromagnetic waves is much less than the proton gyro- 
frequency (for L>5, wave frequency is on the order of or greater than 
proton gyro frequency) and should be mostly governed by hydromagnetic 
instabilities. 
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CHAPTER V 


CONCLUSIONS AND RECOMMENDATIONS 

This study suggests that a theoretical model be used to investigate 
the dynamical characteristics of the wave-particle Interaction in space 
plasma. The formalism of this model is based on the Boltzmann kinetic 
equation. More specifically, the present model dealt only with the 
medium which is in a translational region (i.e., from collisional to 
collisionless regions), namely, the electron equation can be represented 
by fluid equation and ion equation governed by Boltzmann equation. The 
criteria for the validity of this approach is based on the ratio of wave 
frequency and collision frequency. A detailed account of this discussion 
is included in Chapter II. 

This theoretical model has been applied to study the hydromagnet ic 
instabilities in the magnetosphere, in which the instability criteria for 
hydromagnetic wave is established in the transitional region of the 
magnetosphere. 

Possible mechanisms for the firehose Instabilities based on magneto- 
spheric conditions for both quiet and disturbed cases are discussed. It 
is found that the B can be reached to the order of unity in the magneto- 
sphere, then the dissipation rate of hydromagnetic waves jumps to 10“^. 

This gives a magnetospheric heating rate through the damping of hydromagnetic 
waves with amplitude of 10 gamma to be on the order between 10”*^ and 10”^ 
erg-cm" -sec , which is on the same order of EUV solar radiation. Hence, 
we may conclude that the dissipation of hydromagnetic waves can contribute 
to the magnetospheric heating. 
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Further applications of this model will be on the numerical studies 
of various solutions corresponding to those magnetospheric observations j 
and to investigate the region of validity of the present theoretical model. 
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FIGURE CAPTION 


Figure 1 EhAmping rate of hydromagnetic wave for - 1.0 (magneto- 
sphere and Bi = 10“^ (ionosphere) at (jOT^ (the ratio of 
collision time to wave period) =0.5 with <T^> = <Te>. 

This damping rate plotted is based on the numerical results 
of dispersion relation (2-7), 




